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ABSTRACT 


Design  of  experiments  for  estimation  of  parameters  in  non-linear  models 
is  studied  in  a  Bayesian  framework,  with  the  objective  of  maximization  of  the 
anticipated  Fisher  information.  Two  stage  optimal  designs  are  proposed  in 
attribute  life  testing  situations. 

1.  Introduction 

Consider  a  system  of  N  components  working  independently  and  having 
identical  cumulative  distribution  functions  (c.d.f.)  of  the  time  till  failure 
F  C t ;  6 ) .  F  is  a  known  function  and  0  is  an  unknown  parameter  belonging  to  a 
parameter  space  0.  It  is  assumed  that  the  number  of  components,  N,  is  a  fixed 
positive  integer.  The  components  fail  randomly  at  unobservable  times.  We 
inspect  the  system  after  x  units  of  time  and  count  the  number  of  failed  compon¬ 
ents.  The  replacement  of  the  components  could  follow  one  of  the  following  two 
policies : 

(A)  Only  failed  components  are  replaced  at  each  inspection. 

(B)  All  items  in  the  system  are  replaced  at  each  inspection  (frequent 
replacement  policy,  or  block  replacement  policy). 


(*)Part  of  the  Ph.D.  thesis  in  department  of  Mathematics  and  Statistics, 
Case  Western  Reserve  University,  Cleveland,  Ohio.  Supported  by  ONR 
Contract  N00014-80-C-0325  (NR  042-276)  at  Virginia  Polytechnic  Institute 
and  State  University 
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The  policy  to  ascribe  depends  upon  the  type  of  system  under  consideration. 

For  example,  policy  (B)  is  preferred  when  it  costs  more  to  inspect  and  change 
only  failed  components  as  compared  to  changing  the  whole  system.  Such  examples 
are  encountered  in  the  replacement  of  street  light  bulbs,  etc.  Sometimes  it  is 
practically  impossible  to  change  only  the  failed  components  without  effecting 
the  whole  system.  Systems  composed  of  transistors  built  in  modules. 

For  more  applications  see  Barlow  and  Proschan  (1967).  In  quantal  response 
bioassay  studies  policy  (B)  is  followed,  where;  after  experimentation,  the 
whole  batch  of  experimental  units  (mice,  fish,  etc.)  is  replaced  by  a  new  one. 

Finney  (1978)  gives  an  exhaustive  reference  list  of  bioassay  studies  of  this 
kind. 

Let  J(x, ) ,J(x_) , . . . ,J(x  ),...  denote  the  number  of  components  failing 
1  t  n  , 

n-1  n 

during  the  intervals  (0,x  ),  (0,x  ),  (x. ,x  +x_) , . . , (  £  x  ,  £  x.),...  . 

1  1  lit  j.j  i  1 

Intuitively,  we  would  like  to  use  the  information  ( J(x, ) , . . . , J(x  ),x, ,...,x  ) 

Inin 

to  define  x^+i  so  that  J(xn+^)  will  provide  as  much  information  on  0  as  possible. 

To  define  the  best  or  optimal  interinspection  time  at  the  (n+l)st  stage,  we 

shall  use  the  criterion  of  maximizing  the  conditional  Fisher  information  about 

0,  given  (J(x, ) , . . . , J(x  ),x,,...,x  ).  More  specifically  let  F  denote  the  sigma 
i  n  l  n  n 

algebra  generated  bv  (J(x,), . . . ,J(x  ),x,,...,x  )  and  let  l(0;x  ..IF  )  denote  the 
'l  n  l  n  n+x  n 

conditional  Fisher  information  of  0  at  the  (n+l)st  stage  given  F  .  Generally, 
I(0;xn+jjF^)  depends  on  0.  Hence,  the  optimal  value  of  xn+^  is  a  function  of  the 
unknown  parameter  0.  This  problem  can  be  overcome  by  changing  the  criterion  of 
optimality  in  a  suitable  manner. 

In  a  Bayesian  framework  we  assume  that  0  is  a  random  variable  having  a  specified 

distribution  function,  called  the  prior  distribution.  The  Bayesian  interinspec- 
* 

tlon  time  xn+^  i*  a  number  maximizing  the  predictive  Fisher  information  i.e., 

Kll(0;x  ,|F  M  >  EU(0;x|F  )) 
n+1 1  n  -  n 


(1.1) 
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for  n=l,2 , . . . ,  and  V  x  e  X,  where  the  expectation  is  with  respect  to  the  posterior 
distribution  of  0  given  F^  and  X  is  an  appropriate  design  space,  x^  is  defined 
as  the  value  belonging  to  X  such  that 

E{I (0 ;x*) }  >  E{ I (0 ;x) }  V  x  e  X  ,  (1.2) 

where  the  expectation  is  taken  with  respect  to  the  prior  distribution  of  0. 

Zacks  (1973,1977)  discussed  this  problem  when  the  time-till  failure  follows 
an  exponential  distribution 

F(x,0)  =  1  -  exp{-0x} , 

and  0  follows  a  gamma  prior  distribution.  This  case  will  be  studied  along  with 
some  other  distribution  functions  F(x,0)  and  different  prior  distributions  of  0. 

It  is  readily  seen  that  J(x^)  is  a  binomial  random  variable  with  parameters 
N  and  F(x^;0).  Accordingly,  the  Fisher  information  function  of  F(x^;0)  given 
Xj  is 

I (F(x^ ; 0) )  =  N/{F(x1;0)(l-F(x1;0))}.  (1.3) 

Therefore  the  Fisher  information  of  0  given  x^  is  (Khan  (1 980)) 

KfijXj)  =  I(F(*1;e))(^F(x1ie))2  .  (1.4) 

Unfortunately,  the  conditional  distribution  of  J(x^)  given  F,  is  not  necessarily 

binomial  under  replacement  policy  A,  unless  F(x;0)  is  a  negative  exponential.  This 

complication  arises  due  to  the  fact  that  at  the  2nd  stage  we  have  two  kinds  of 

failure  distributions  -  for  those  components  which  failed  in  the  previous 

interval  and  were  replaced  by  identical  components  we  still  have  the  failure 

distribution  F(x.,;0),  however,  for  those  components  which  did  not  fail  during 
* 

(0,Xj)  the  failure  distribution  Is 
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G(x2;x*,0)  =  [F(x*  +  x2;0)  -  F(x*;0)]/[1  -  F(x*;6)]  . 

* 

Clearly,  if  F  is  negative  exponential  C(x2;x^;0)  =  F(x2;0).  So 

j(x2)|F1  =  Y  +  Z 

where  Y  ~  B(J(x*),  F(x2;0)) 

Z  ~  B(N  -  J(x1),  G(x2;x1,0)) 

and  Y,  Z  are  condi tionallly  independent. 

Now  I(0;x2|F^)  does  not  have  the  same  form  as  (1.4)  for  general  failure  dis¬ 
tribution  F(x;0).  This  problem,  however,  can  be  overcome  by  the  method  used 
by  Zacks  and  Fenske  (1973).  If  F(x,0)  is  the  negative  exponential  distribution, 
then 

I(0,x2jF1)  =  1(0, x2) 

* 

and  the  second  stage  optimal  interinspection  time  x2  is  the  value  maximizing 

E{I(0;x2)} 

where  the  expectation  is  with  respect  to  the  posterior  distribution  of  0  given 
F  .  Under  replacement  policy  B,  on  the  other  hand, 

J ( x  )|F  ~  B(N,F(x  ;0)  ;  n=l,2,...  , 

n  n- 1  n 

for  all  distributions  F(x;0).  Hence, 

I(9;xJFn-l}  =  I<F(V0))<80  F(xn’e))2  *  (1.5) 

Thus,  under  replacement  policy  B  the  problem  has  the  same  structure  in  all  the 
stages.  In  the  case  of  negative  exponetial  distribution  the  problem  remains 
the  same  in  both  kind  of  replacement  policies. 
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2.  One-Stage  Designs 

In  this  section  we  shall  discuss  the  determination  of  interinspection  times 
for  the  following  cumulative  distribution  functions: 

(i)  F(x,0)  =  1  -  exp{-x/0};  the  negative  exponential  distribution. 

(ii)  F(x,0)  =  2(1  +  exp{-0x})  *  -  1  ;  x  >  0  ,  the  truncated  logistic  distribution. 

(iii)  F(x;0)  =  (e-1)  (e  •  exp{-e  1-1}  ;  x  >  0  ;  the  truncated  extreme 
value  distribution. 

(iv)  F(x,0)  having  a  symmetric  density  function  and  0  is  the  shift  parameter. 

These  examples  will  suffice  to  show  the  complexity  of  the  algebraic  manipulations 
involved  for  the  solution  of  this  problem.  Nevertheless  the  method  is  straight¬ 
forward  and  could  be  applied  to  any  F(x;0). 

(i)  The  negative  exponential  distribution: 

If  F (x, 0)  =  1  -  exp{-x/0};  x  >  0,  0  >  0,  then  by  equation  (1.4) 

l(0;x)  =  Nx2/{0**(exp(x/0)-l ) }  .  (2.1) 

The  design  level  x  maximizing  (2.1)  is  the  solution  of  the  equation 

exp(x/0){2  -  x/0)  -2=0, 

which  is,  x^  =  1.5936  0. 

The  prior  information  with  respect  to  a  prior  distribution  G(0)  is  defined 
as  the  prior  expectation  of  l(0;x)  and  is  given  by 

E{  1(0  ;x) }  =  Nx2E{e"x/O/{04(l-e-x/e)H 

cu 

=  Nx2E(  I  exp{-kx/0}/0^} 
i=l 

GO 

=  Nx2  Z  E{exp{-kx/0 } /Q* ) 
k=*l 


(2.2) 
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The  exchange  of  summation  and  expectation  in  (2.2)  is  permissible,  since  the 
function  is  non-negative. 

The  Bayes  optimal  design  level,  x^(G),  is  a  value  of  x  satisfying 
“  M  -kx/G 

E  r  - — 7 —  O  -  kx/20)  dC(0)  =  0  (2.3) 

k=i  o  e 

provided  the  left  hand  side  of  (2.3)  converges  uniformly  in  x. 

Consider  the  case  when  G(0)  is  the  inverted  gamma  distribution  with  para¬ 
meters  (A,m).  The  corresponding  density  function  g(0)  is  given  by 

A*  -m-5 

g(G)  =  ^  P  exp(-X/0)  ;  0  >  0  . 

For  this  prior  distribution  (2.2)  reduced  to 

cc 

E(I(0;x)  }  =  cx^  T.  (X  +  xk)  m  ^ 
k=l 

oo 

=  c'x2  T.  (1  +  xk/X)_m~A 
k=l 

when  c  and  c'  are  appropriate  constants. 

Without  loss  of  generality,  we  discuss  the  problem  of  choosing  an  x  to 
maximize 

00 

f(x)  =  Z  x2/(l  +  xk)11^^ 
k=l 

* 

The  solution  Xq  will  provide  the  optimal  design  level  x^  for  all  X  by  the 
relationship 

x*  =  XxQ  . 

Let , 


g(k,x)  =  x2/(l  +  xk)"^ 
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Lemma  2.1. 

00 

(i)  I  g(k,x)  is  uniformly  convergent  in  X  e  (0,®). 
k=l 


(ii)  ^  f(x)  =  Z  ^  g(k,x)  . 
k=l 


(iii)  Xq  ~  1.5936/(m+2)  for  large  values  of  m. 
Proof. 

Each  function  g(k,x)  attains  a  unique  maximum  at 


x^  =  2/(k(mf2))  ,  k=l,2,...  . 


Let, 


^  =  gCk^)  =  sup  g(k,x)  ,  k=l  ,2 , . 


It  follows  that 


l  M.  ~  Z 
k=l  *  k=l  k2 


Thus  (i)  is  proved  by  the  Weirstrass  M-test  (Widder  1961). 
(ii)  follows. 

To  prove  (iii)  we  notice  that 

g(k,x)  >  g(k+l,x)  V  x  e  (0,®),  k=l , 2 , . . . 

Since  g(k,x)  attains  its  maximum  value  at  x,^  =  2/(nrf2)k,  it 
for  each  k=l,2, .  . . 

r 


By  similar  arguments 


readily  follows  that 


g(k,x)  = 


+  for  x  e  [2/(m+2)k  ,  ®) 
f  for  x  e  (0,2/ (nrt-2)k) 


(2.4) 
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Indeed , 


and 


x  e  (0,  2/ (nri-2) ) 
m,  u 

g(k,x)  £  g(l,x)  ;  k=l  ,2, , 
\  1  Xx  »  k=l  ,2, , 


According  to  (ii)  one  can  differentiate  under  the  summation  and 


—  f(x)  =0  is  equivalent  to 


£  {2  -  (m+2)kx)  (1  +  kx)  =  0 

k=J 


(2.5) 


Let  x  =  a  (nH-2)  ^  where  for  each  m,  0  <  a  <2.  The  values  of  x  (or  a  )  can 
mm  m  mm 

be  determined  numerically,  for  each  m  from  equation  (2.5).  We  show  now  that  f.r 
large  values  of  m  the  solution  has  a  simple  approximation.  For  this  purpose  v 
establish  first  that  a  is  a  convergent  sequence.  Indeed,  there  exists  a  sub¬ 
sequence  and  a  limit  point  a  such  that 


I 

I 


Hence , 


a  •+  a  e  (0,2]  as  v  +«> 
m 

v 
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Taking  the  limit  as  v  -*■  +»,  one  obtains  the  equation 

2  -  a  =  E  (ak-2)e"a(k_1)  , 
k=2 

or 


2  -  a  =  e 


-a 


a(2-e~a)  2 


l/l  *~n .  2  i  u 

L (1-e  )  1  -  e 


this  equation  is  further  reduced  to 


2-a-2e=0. 


i 


'  i 

)  i 

V, 

!; 


1 


ei 


The  solution  of  this  equation  is  approximately  a  =  1.5936.  i.e.,  a  has  an 

m 

unique  limit  point.  Hence,  x  .  ~  1 . 5936/ (m+2) ,  as  m-«°. 

m,  U 

(Q.E.D.) 

•k 

Some  values  of  a  ,  x  n  and  x  for  various  values  of  m  are  given  in  Table 
m  m,U  m 


The  results  obtained  so  far  are  summarized  in  the  following  theorem. 


Table  2.1 


Values  of  a  ,  x  and  x  „  for  A=1 
mm  m,  0 


m 

a 

m 

* 

X 

tn 

Xm,  0 

0 

1.530 

.765 

.797 

1 

1.549 

.516 

.531 

2 

1.556 

.389 

.398 

3 

1.560 

.318 

.318 

4 

1.564 

.261 

.266 

5 

1.567 

.224 

.228 

6 

1.568 

.196 

.199 

7 

1.573 

.175 

.177 

8 

1.577 

.158 

.159 

9 

1.580 

.144 

.145 

15 

1.582 

.093 

.094 

20 

1.584 

.072 

.072 

+oo 

1.594 

0 

0 

1 
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Theorem  2.1. 

If  N  components  are  working  independently  with  identical  failure  distribution 

F(x,0)  =  1  -  exp(-x/0) 

W.e  rc  0  has  the  inverted  gamma  (A,m).  X,m  >  0,  prior  d'str'bntion ,  'hen  the 

» .  .  .  * 

first  stage  optimal  intcrmspeetion  time  x^  ^  is  given  by 

*  - 1 
x  =  A  a  ( m+2 ) 
m,  1  m 

where  the  values  of  a  are  given  in  Table  2.1. 

m 

v.  x  1 )  I truncated  logistic  distribution  . 

We  consider  the  logistic  distribution  truncated  to  the  left,  i.e., 

F(x, 6)  =  2(1  +  exp(-Qx) )_1  -  1  ;  x  >  0  ,  0  >  0  . 

In  this  case  the  failure  (hazard)  rate  function, 

h(x)  =  0(1  +  exp(-Ox))  1  , 

is  an  increasing  function  of  x. 

Also  according  to  (1.4)  the  Fisher  information  function  of  0  given  x  is 


1(6, x)  = 


_ 2Nx  exp(-Ox) _ 

l-exp(-Ox) } 1 l+exp(-0x) }i 


The  design  level  x  maximizing  (2.6)  is  the  unique  solution  of  the  equation 

o  q  a  / o  “9x.  /  20X  . -1 

2  -  0x  =  0x(3  -  e  ) (e  -1) 

which  is  x°(ti)  =  2.I6ji/0. 


For  a  given  prior  distribution  of  0,  the  expected  Fisher  information  function 


is  given  by 
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E{I(6,x)}  =  2Nx  E  { 


,  -0x.  ,  -6x.  2 

(1-e  ) (1+e  ) 


CO  oo 

=  2NX2  1  E  (-l)j(j+l)E{e-(k+j)x9} 
k=l  j=0 


(2.7) 


=  2Nx2  E  E  (-l)J(j+l)M  (-(k+j)x) 
k=l  j=0 


where  M  denotes  the  Laplace  Steiltjes  transform  of  the  c.d.f.  of  0.  It  is 

o 

rt 

difficult  to  characterize  the  design  level  x  which  maximizes  (2.7)  in  the 
general  case.  Therefore,  consider  the  special  case  where  0  has  a  prior  gamma 
distribution,  G(A,m),  with  m  >  2. 

In  this  case 

M0(-(k+j)x)  =  (1  +  x(k+j)/\)  m  . 

Hence,  without  loss  of  generality  we  can  assume  that  \  =  1,  and  consider,  for 
m  >  2,  the  expression 

CO  oo 

E{  I  (6 ,  x)  =  CX2  E  E  (-l)j(j+l)U  +  x(k+j))"m  (2.8) 

k=l  j=0 

Differentiating  with  respect  to  x  under  the  summation  signs  we  obtain 


p-  E{  1(0 ;  x)  }  =  c  “  I  (-l)j(j+l)  . 

k=l  j=0  (l+x(k+j))m  1 


Let  a  =  x  (m-2)  where  x  is  the  solution  of  the  equation 
m  m  m 


2  -  a  (k+j) 
m 


E  E  (-l)J(j+l) - - - t-  ,-r-  =  0  V  m  , 

k=l  j=l  (1  +  (k+j) a  (m-2)  ) 

m 


(2.9) 


obtained  by  equating  the  partial  derivative  to  zero.  To  see  why  is  a  con¬ 
vergent  sequence  as  we  note  that  for  large  m 
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(1  +  (k+j)u  (m-2)  •^)nrt’^ 
m 


a  (k+j) 
m  J 
e 


i.e.  r. 

k=l 


00  .  2  -  a  (k+j ) 

T.  (-l)J(j+l)  — - -- - 

.1=0  ^  m(k+j) 


0  for  large  m 


=>2(1 


-2a  -a  -3a  -2a 

/ 1  I  rn «.  / 1  i  o  ni  .  ni» 

e  )(1  +  e  )  -  a  (1  +  2e  +  e  ) 

m 


0  for  large  m 


=>  a  ~  2.1651  for  large  m 
m 

Hence,  a  -*■  2.1651  as  m -*■+“’.  As  before,  when  X  ^  1 ,  we  have 
m 

x  =  Xa  (m-2) 
m  m 

In  the  following  table  we  give  the  optimal  solution  in  terms  of  a  and  x  for 

m  m 

small  values  of  m,  as  determined  by  numerical  solution  of  (2.9). 


Table  2.2 
X  =  1 


m 

a 

m 

* 

X 

m 

x  =2 . 1651 (m-2)-1 
m 

3 

2.04687 

2.04687 

2.1651 

4 

2.06844 

1 .03422 

1.0825 

5 

2.08836 

0.69612 

0.7217 

b 

2.10172 

0.52543 

0.5413 

9 

2.12331 

0.30333 

0.3093 

10 

2.12768 

0.26596 

0.2706 

15 

2.14032 

0.16464 

0.1665 

+*> 

2.16509 

0.0 

0.0 

The 

truncated  extreme  value 

distribution . 

The  truncated  extreme  value  distribution  is  given  by 

ir '  /. . ') )  -  (o-l)  '  (e  exp(o  ^ X)  -  1 )  ;  x  ■  0  ,  o  0  . 


1 

J 
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The  Fisher  information  function  of  0  at  a  design  level  x  is 

K6, x)  =  ^_xi_exp(-2e-6x)_e_xR(,-26_x^_  .  (2.10) 

[e  exp(e  x)-l ] [l-exp(-e  °X)  ] 

The  design  level  x  maximizing  (2.10)  satisfies  the  equation 

2(1  +  l/f»tz)/z  =  [(l+e)exp(-z)-2]/[(l-e  2)(e  e  Z-l)  ]  , 

where  z  =  exp(-Ox).  The  solution  of  this  equation  yields 

x  =  -Cn(. 144188) /0  =  1.9366/0  .  (2.11) 

In  analogy  to  the  previous  two  cases  we  may  conjecture  that  if  the  orior 
distribution  of  0  is  gamma  G(A,m),  the  optimal  design  level,  at  which  E{I(0;x)} 
is  maximized,  is 

x^  ~  A(1.9366)/(m-2)  for  large  m. 


Theorem  2.2 

Suppose  that  the  Fisher  information  function  is  of  the  form  1(0 ;x)  = 

I  a^xPexp{-6kx} ;  p  j  0,  where  a^  do  not  depend  on  0  and  x;  and  that  I(6;x)  is 

k  .  .  0 

maximized  at  x  (0)  =  a/0.  Furthermore ,  assume  that 


(i)  Z  a  x^  ^(p-xk0)e  is  uniformly  convergent  in  x, 

k  k 

, .  . ,  _  Am-1  -(X+xk)9  ..  .  . 

(n)  Z  a,  0  e  is  uniformly  convergent  in  0, 

k  k 


r  p-1  (p-(m-p)xk)  ..  ,  _  . 

(m)  tax  --*■ - p  -  is  uniformly  convergent  in  x, 

k  k  (l+xk)11”"1 


then  if  the  prior  distribution  of  0  is  the  gamma  G(A,m),  E{I(0;x)}  is  maximized 
at  Aa/(m-p)  for  large  values  of  m. 


M 


Proof.  Since  l(0;x)  is  maximized  at  x  =  a/0,  therefore 


9x 


l(0;x)  =  0  for  x  =  a/0  . 


(2.12) 


—  0 

It  follows  from  (i)  that  E  a  (p-xkB)e  =  0  for  x  -  a/0 

k  k 

assumption  (ii)  we  obtain  that 


From 


E{I(Q;x)}  =  E  a  xP  Ele  I 
k  k 


E  a  xP(l  +  xk/X)  m  . 
k  K 


(2.13) 


tlv.'-it  loss  of  generality,  assume  that  X  =  1. 
Furthermore,  from  (iii)  we  get 


ri  E{I(€; x) )  -  x^'V  P  -  ■ 

’*  k  k  (l+xk)”*1 


(2.14) 


net  x  be  the  value  of  x  for  which  (2.12)  is  equal  to  zero,  and  let  a  =  (m-p)x 
m  m  m 


We  thus  obtain  the  equation 


Ea  (p-a  .  )(1  +  k)  m  1  =  0  V  m. 

,  k  1  mk  m-p 


(2.13) 


Moreover, 


0  =  Ea  (p-a  k) (1  +  — ^— )  m  1 
k  k  m  m-p 


-a  k 

Ea,  (pa  k)e  for  large  m 

,  k  m 
k 


=>  a  ^  a  for  large  m  , 
m 

a  k  a^k 

since  (l  +  — — )  ~  e  .  Hence,  x  ~  Xa/(m-p),  for  large  m. 

m-p  m 


(Q.E.D. 
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Remarks.  1.  For  the  truncated  extreme  value  distribution 


1(6 ;x) 


[e 


„  2  ,  _  -0 x.  -26x 

Nex  exp(-2e  )e _ 

exp(-e  -  1 ] [ 1  -  exp(-e  ®X) ] 


could  be  written  as 


I (6 ; x)  =  Nx 


j,k,£,q,v=0 


aj,k,i,a,ve 


-(2+j+Jl+v)0x 


where 


=  (-2)j(-k)a(q+l)Vq 

j  ,  k ,  9  ,  q ,  v  j  !  9 !  q!  v ! 

2.  The  conditions  (i),  (ii),  (iii),  of  Theorem  2.2  could  be  relaxed  after 

. i  ;  ■ .  !  i  .  :  divergence  of  the  series 


la  exp(-Bxk)  . 
k  k 

The  following  theorem  generalizes  the  last  theorem  for  the  prior  distribu¬ 
tions  belonging  to  the  class  of  infinitely  divisible  distributions. 


Theorem  2.3. 

If  the  Fisher  information  function  of  6  is  of  the  following  form 

l(0;x)  =  la  xPexp(-8xk)  ;  p  0  , 
k  * 

where  ak  do  not  depend  on  8  or  x;  and  if  I(6;x)  is  maximized  at  the  design  level 
a/6  and  6  has  a  prior  distribution  belonging  to  the  class  of  infinitely  divisible 
distributions  with  parameter  m,  then  E{I(0;x)}  is  maximized  at  the  design  level 
x^  ~  a/E{6)  for  large  m  provided 

(i)  Ta  x17-1  (p-xkO)exp(-0xk)  is  uniformly  convergent  in  x. 
k  k 

(ii)  5:akg(0)  exp(-Oxk)  is  uniformly  convergent  in  0. 

(iii)  l^kX^- *E( (p-xOk)exp(-Bxk) }  is  uniformly  convergent  in  x. 
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(iv)  m  x  ia  a  bounded  sequence 

m 


Proof . 


Since  1(6 ;x)  is  maximized  at  x  =  a/0, 


Now,  according  to  (ii). 


-  l(i.!;>)  =  0  for  x  =  u/0 

d  X 


E{  I  (0 ;  x)  1  *  la,  xpE«e  l‘Xk( 
k  k 


(2 . 16) 


=  Ea  x%(xk)  , 
k  k 


where  <j)(t)  =  exp  (— n>4»  (t ) ) ,  such  that  tj>(0)  0. 


Hence,  by  (iii) 


E{I(6;x  )}  =  La  {p  +  mxk  ip  (xk)  !  exp  (-m<|i(xk)  }  =  0. 
3x  m  ,  k 

k 

Let  a  =  m  x  .  We  thus  have  the  equation 
mm 


La.  {p  +  ka  \p  ( —  ) }exp(-m  ip ( - ))  =0  V  m  , 

,  k  mm  m 

k 


(2.17) 


since  a  =  mx  is  bounded,  there  exists  a  subsequence  m  such  that  m  •  x 
mm  v  v  m 

v 

converges  to  a  finite  limit,  say  a',  as  v  ►  +®.  Hence,  x  +  0  as  v-H-®. 

m 


Expanding, 


we  have 


p(x)  =  i|» (0)  +  x^'(0)  +  — j  ♦"(0)  +  ... 


*(— )  =  —  k  **  (0)  +  (— )2  k4  *"(0)  +  ... 


r'  m  '  m  T  ' '  '  m'  2 

a  (a  k)^ 

-iw|>( —  k)  =  -a  k  if)’(O) - r - 

ni  m  Z  m 


*"(0)  +  ...  . 


^ - 


1 
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Therefore , 


Or, 


k  a 


m 

exp(-m  iji( - ))  -*■  exp(-a'k0' (0))  as  v  *  +»  . 

v  m 

v 


Sa  Ip  -  a'k(-(|,' (0))  lexp(-«'M''(0))  =  0 
k 


Hence,  o'  =  (0))  and  a  •  a/ i -tji '  (0)  I ,  which  implies  that 

m 

x  ~  in  /  E  {  0  1  for  largo  values  of  m  , 
m 


(Q.E.D.) 


Remarks : 

(i)  As  the  previous  examples  show,  the  approximation  of  x  is  better  if 

m 

-  -mil  'fOl  i  s  r^rlaoed  hv  -  (m-p)i)  ' (0) . 

(ii)  Zacks  (1973,1977)  considered  the  case  when 

F(x;0)  =  1  -  exp(-Bx) 
wi th  D  ~  C( A ,m)  ;  m  •  2  , 
by  numerical  results  it  was  conjectured  that 


x  =  1.5° 36  \/(m-2)  . 
m 


By  the  above  theorem  we  see  that  this  is  an  aysmptotic  result  for  large 

values  of  m.  However  for  small  values  of  m,  x  =  r>  A/(m-2)  where  a  are  given  in 

mm  m 

Table  2.1. 

( i v )  F(x,0)  having  symme trie  density  function . 

So  far  we  have  discussed  the  cases  when  the  parameter  of  interest  is  a 
scale  parameter.  However,  some  times  we  are  interested  in  the  shift  parameter. 

Consider  the  case  where  F(x;8)  is  the  logistic  distribution  with  a  shift 
parameter 

F(x;0)  =  (1  +  exp(-x+0)J  ^  ,  -*»  <  x  •  ",  *-  0  <  »  . 
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TTu  '‘flirt;  1(0;  x)  is  maximized  at  the  median  of  F(x;G),  furthermore,  since  F(x;0) 
is  symmetric,  x^(G)  =  0.  The  Bayesian  estimation  of  the  median  of  the  logistic 
distribution  was  discussed  by  Freeman  (1970).  However,  dvnamic  programming 
methods  were  used  to  obtain  sequentially  optimal  designs,  up  to  three  stages, 
when  the  prior  distribution  of  0  was  taken  to  be  the  conjugate  prior.  We  shall 
consider  this  problem  from  the  point  of  view  of  maximizing  the  Fisher  information 
For  any  symmetric  c.d.f.  F(x;0),  the  Fisher  information  function  of  0  given 

x  i  s 

l(0;x)  «  N{^|  F(x;9)}2/[F(x;0)(1  -  F(x;0))] 

=  N  f ^(x; 0) / [F(x; 0) ( 1  -  F(x;0))]  ,  (2.18) 


wlie re  f(x;0)  is  the  density  of  F(x;9). 

Consider  the  case  when  f(x;0)  is  symmetric  about  x  =  0.  Note  that  f(x;0) 
is  defined  on  the  real  line  R.  Otherwise,  the  Fisher  information  function  of  0 
does  not  exist.  Furthermore, 


—  f(x;fi)  ->-0  as 
dx 


We  also  note  that  in  order  that  I(o;x)  be  maximized  at  x  =  0  we  need  the  necessary 
conditions  that 

(i)  f ' (fa; e)  =  0 
(ii)  f"(e;e)  <  o 


where  prime  denotes  differentiation  with  respect  to  x. 

This  implies  that  the  density  f(x;0)  is  maximized  at  x  =  0.  Moreover,  if 
f(x:o)  attains  a  minimum  value  at  x  =  o,  then  I(o;x)  is  minimized  in  a  neighbor¬ 


hood  around  0. 


£ 
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Accordingly,  we  consider  only  those  symmetric  densities  which  are  bell 
shaped  (e.g.,  logistic,  normal,  etc.,  commonly  known  as  Logit  and  Probit  models 
in  bioassay). 


Theorem  2.4. 

Let  1(0 ;x)  be  the  Fisher  information  function  of  the  shift  parameter  6  at 
the  design  level  x  for  distributions  F(x;0)  symmctri.e  about  0.  If 


(i)  f '  ( x ; 0 )  £  0  V  x  >_  0, 

(ii)  f"(0,0)  +  2  f(3)(0,0)  <  0  . 


Then  x  =  0  in  a  point  of  maxi-ma  of  I  ( n ;  x) . 

r  •  ?  ,,  •  r 


(iii)  /  f(t,0)dt  <  1/2  [  1  -  (f(x,O)/f(0,0))2]1/2  V  x  '■  6 . 
0 


Then 


Proof . 


Since 


we  obtain 


l(0;x)  attains  its  maximum  value  at.  x  =  0. 


I(0;x)  =  N  f  ( x; 0) / | F(x; 0) ( 1  -  F(x;0))]  . 


1(0; x)  ~  N  f  (x;0)/{]  -  F(x;0)]  as  x  -*•+<*>  , 


lira  I (0 ;x)  =  -2  N  lim  f'(x;0)  =  0  . 
xH-“ 


Similarly  lim  l(0;x)  =  0.  Thus,  there  exists  a  design  level  x  such  that  l(0;x) 

X+-«i> 

is  maximi zed .  Now  differentiating  logl(0;x)  and  equating  to  zero  we  get 


2f,(x,0)  f (x;0)  f(x;9) 

f (x; 0)  '  F(x;0)  1  -  F(x;0) 
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If  x*  0,  then  £'(x;0)  -  0  and  F(x;0)  =  1  -  F(x;0).  Therefore  —  log  K9;x)|x»e 
=  0  1 . e . ,  x  =  0  is  a  point  of  extremum  of  I(9;x). 

By  twice  differentiating  the  log  l(6;x)  we  get 


i 

VI 


■ 

l 


log  1(6 ;x)  »  2 


3x 


f (x:6)f"(x;6)  -  (f *(x;6))' 
f2(x;0) 


(1-F(x:0))f ’ (x;6)  +  f2(x;0)  F(x:6)f *(x;6)  +  f  (x;6) 


[\  -  F(x;9) 1 


(F(x;0)) 


Now, 


we  ge  t 


f '  (6  ; 0)  -0,1-  F(6, 0)  =  F(0;O)  =  1/2 


log  l(0;x) 


3x 


_  -  2  U  f2(e’e>  <  0  ^  (ii)  ’ 
x-0  f(0,0) 


.e.,  I(0;x)  attains  a  local  maximum  value  at  x  =  0. 

Now  by  (ill) ,  we  get 

x  2 

F(x;0) (1  -  F(x;8)>  -  (h  -  (J  f(t;0)dt)  ) 

0 

>  v  x  >  0 


therefore, 


Jr, 


f  (6,0) 


1(0,9)  >  1(0, x)  V  x  >  9. 

Siine  f 2 ( x, 0 )  is  symmetric  about  x  =  0  one  can  imply  that  1(9, x)  is  also 
symmetric  about  x  =  0.  Hence, 

1(0,0)  >  1(0, x)  V  x  i  0  . 


Therefore  l(9,x)  attains  its  maximum  at  x  -  6, 


(Q.E.D.) 
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*  ;i  +  b 
1 .  e .  x  =  — ^ —  • 

Hence  the  first  stage  optimal  Bayes  design  level  is  the  median  of  the 
uniform  distribution.  We  can  generalize  the  above  result  to  any  symmetric 
prior  distribution  as  follows. 

In  general  let  g(t)  be  a  symmetric  bell  shaped  density  symmetric  about  6. 
Then  for  the  prior  distribution  of  8  having  density  g(6-X)  defined  over  the 
real  line,  the  expected  Fisher  information  function  of  8  is 

E ( I (9 . x) )  -  »  /  ft --yi  ECS-X) 


d0 


CO  2 

■  "  l  FCjfo'-HfyTT  dy 


=  /  Ky)  g(x-y-X)  dy 


where  I(y)  =  f  (y)/F(y){l  -  F(y)} 


For  the  sake  of  maximization  of  E{I(0,x)}  we  see  that  without  loss  of  generality 

we  can  assume  that  X  =  0.  For  the  following  theorem  we  assume  that  F(x,0)  = 
x 

/  f ( t  —  *"> )  dt  when  f  satisfies  the  conditions  of  the  previous  theorem  such  that 

—  n 

I(fl,x)  is  also  bell  shaped. 


Theorem  2,5. 

Li  t  I (0 , x)  he  the  Fisher  information  function  of  0.  Suppose  that  0  has 

a  prior  distribution  G(0-X)  defined  on  R,  such  that  its  p.d.f.  g  is  symmetric 

about  the  median  X.  Then  x  =  X  is  the  point  of  maxima  of  E{I(0,x)}  provided 

00 

/  g"(e)  d0  <_  o  .  (2.19) 
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Proof.  Assume  that  X  =  0.  Then 


more  tore. 


CO 

E{I(9 ;x) }  =  /  I (y)  g(x-y)  dy  . 

—CO 


E{I(B;x)  =  /  I(y)  gl'x~y)  dy 


=  0  if  x  =  0  =  }  . 


This  follows  by  considering  the  fact  that  g'(-)  is  an  odd  function  and  l(y) 
is  symmetric  about  y  =  0.  Hence,  x  =  X  =  0  is  an  extremal  point  of  the  expected 
/.Isa i  •'  j.ji  orm-itj-on  function.  FurtUennoie, 

a2 

E{  I  (0  ;  x)  }  =  /  I(y)  g"(x-y)  dy  . 

3x  -<» 

Lev  -HC  denote  the  inflexion  points  of  the  density  gf  •  )  •  Hence, 

2  "» 

E{  1  (d ,  x) }  l  =  /  I  (y)  g"(-y)  dy 
*v“  x-n 

00 

=  /  I(y)  g"(y1  dy  because  g"  is  an  even  function. 

—  oo 

-K  K 

-  (/  +  f  +  j  X  Tfv)  g"(v>  dy 

-K  K 

“  K  K 

=2  /  I (y)  g"(y)  dy  +  /  l(y)  g"(y)  dy 

-co  -K 

-K  K 

<  2  l(K)  f  g" (y )  dy  +  T (K)  /  g"(y)  dy  . 

~°°  -K 


I'h i  -  follows  by  noting  that 


g"<e) 


<  0  V  0  <  (-K.K) 
0  V  0  /  (-K,K) 


a2 

Hence,  — r  Ell(0,x)}  |  <  I(K)  /  g"(y)  dy  0  =  x*  =  A  is  the  point  of  maxima 

3x  x=0 

the  expected  Fisher  information  function. 

(Q. 

Example 

If  0  ~  N(0,cr~) ,  then 

at 

/  g"(y)  dy  =  E(e2/o‘!'  -  ~  ) 

-o>  O 

=  0 

Therefore  the  normal  density  satisfies  the  condition  of  Lhe  above  thuwi. - 


3.  Optimal  Design  of  Two  Consecutive  Design  Levels 

W  sh:  11  ,.cn  iaor  the  cv-<v.c,,f ;  f:i’  •  1  ,  : 

section.  Therefore  the  results  obtained  ate  applicable  under  both  policies 


(A;  a  t;!  1.H)  . 

*  .  * 
If  J  (xj),  J2 (x?)  denote  the  number  rf  components  failing  during  fO,Xj ) 

*  * 

(x^,x2  +  x2)  respectively,  -hen 

J2 (X2> I J 1 ^ x l > ’ xo  ~  H(N,F(x2,6)) 

* 

and  the  conditional  Fisher  information  of  0  given  J^Cx^)  and  *2  is 


1(0, x2) 


N  x„ 


4 

0  (e  “  -  1) 


THus,  if  9  has  an  inverted  gamma  (A,m)  prior  distribution,  then  the 

"k 

posterior  distribution  of  6  given  J,  (x^)  =  i  is 

f(0jj)  =  C(x  ,.j)b(j  ,N,F(x  ,0))  tt  e 


L‘  t 


£-D.) 


aii 


,  0  <  e 


25 


It  follows  that 


-^,/0  j  :<1(N-j)^+x2) 

KU(6;x9)  |  J  (x*)  =  j,xj  =  c'(x*  j)x?  f  (1~6  )  GXp(~  ~  6~  '1 

o  e  -x  /e 

(1-0  ) 


dO 


GO  -j 

c'(x  j)x~  1  j:  (J)(-.L)C 
k=l  1=0 


(3.1) 


(A  +  x  (N=j+{,)  +  kx,) 


nH-4 


In  the  special  case  of  j  =  0,  the  optimal  second  stage  interinspection  time, 
* 

say,  q  is  the  point  at  which 


jfc  nv-f-A 

k*l  (\  +  Nx  +  x2k)^ 


is  maximized 


Applying  Lemma  2.1  we  obtain. 


2,0 


(A  +  N  x,  )a 
_  1  m 

n;  +  2 


(3.2) 


where  a  is  given  in  Table  2.1  and  a  1.5936  as  m  -*•  +«’. 
m  m 

* 

Let  x?  .  denote  the  second  stage  optimal  interinspection  time  given  that 

^  9  j 

j  number  of  components  failed  during  (0,x  ).  We  shall  call  an  inspection 
"redundant"  if  the  number  of  components  failed  is  zero. 

k 


Intuitively  x^  ^ 


This  is  because  no  component  failed  during  (0,x^).  So  we  would  give  the  system 
* 

more  than  x^  units  of  time  so  that  some  of  the  components  may  fail  and  avoid 
any  redundant  i:  peer u a:  . 

■k 

Also  if  j(>0)  number  of  components  failed  during  (0,.\j),  then  the  inspection 

* 

was  not  redundant.  So  the  next  interinspection  time  .  should  not  be  as  large 
as  had  I  he  first  inspection  produced  a  redundant  inspection  i.e., 

x9  i  1  x9  n  V  J  =  0 , 1 , 2 , . .  . ,  N  . 


^ - 
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On  the  other  hand,  if  j  =  N,  i.e.,  if  all  the  components  had  failed  at 
the  time  of  inspection,  we  let  the  system  run  too  long  without  any  inspection. 
Therefore, 

*  * 

X2,N  -  X1  ' 

And  if  j(<N)  number  of  components  failed  during  (O.x^),  then  the  2nd  stage 

*  * 

interinspection  time  x„  .  should  not  be  as  small  as  x.  ,,,  i.e., 

2,2  2,N 

X2,N  -  x2,j  V  j  0,1,2 . N* 

We  can  combine  the  above  intuitive  results  as 


y.,  ,,  <  x.  ,  < 
2, A  —  2,N-1  — 


*  *  * 

.  <  X  <  X  <  x0  .  < 

-  2,J0  -  1  -  2 , j rt+l  - 


'0 


-  X2,l  -  X2,0  ’ 


0.3) 


where  is  some  integer  belonging  to  the  set  { 1 , 2 , . . . ,N-1 } .  Again,  intuitively, 
we  feel  that  =  .80(N),  because  the  maximal  Fisher  information  of  6  is  obtained 
at  the  80th  percentile  of  the  exponential  distribution. 

The  following  lemma  verifies  the  above  intuitive  results. 

Lemma  3.1. 


E{I(0,x  |F  )}  is  maximized  at  x  .  where 

^  1  »  J 


2.3 


m+2  +  (N"1)xl 


9  aJX 

jx,/(e 


-  1)} 


(3.4) 


and  a  -*■  1.5936  as  m -»•+'».  (a  are  given  in  Table  2.1). 
m  m 


Proof . 


E{I(8;x  \T  )}  =  C(x  ,j)x  E  I  (,)  (-1) 

k=l  Jl=0 


(X+x*(N-j+Jl)+kx2)nr+'4 
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1  * 

- 7 t  >  where  a„  «  X  +  x, (N  -  j  +  %) . 

nrK»  S,,m  1  J 


*JO  j 

Let  f (x)  =  x2  l  l  dX-l)4  - 

k=l  =0  (a.  +  xk) 

*  > m 
Notice  that  may  depend  on  m.  Approximate  the  function  f(x)  by  using  integrals 

00 

instead  of  the  summation  I  . 

k=l 

Note  that  there  exists  a  constant  c  .  .  such  that 


f(x)  =  x*  £  (j)(-l)  / 


dk 


i.=0V  c  \ '(a  +  xk)®1'4 

m,j,i  l,  m 


Furthermore , 


f(x)  «  x2  l  (h(- n*  ! 


dk 


i=0  c  (a.  +  xk) 

m  i  ,m 


where  c  e  [0,1) •  (We  are  approximating  f(x)  by  choosing  c  independent  of  l 
m  m 

and  j . )  Now, 


dk 


c  (a  +xk)m+4  x(m+3)(a  +  c  x)m+^ 

m  H ,  m  C.  ,m  m 


Therefore, 


f (x)  =  I  (j)(-l)  - - - rrr  . 

i=0  *■  (m+3)  (a  .  +  c  x) 

!t,m  m 


Moreover,  —  f(x)  =  0  implies  the  equation 


j  a  +  c  x  -  c  x(nrl-3) 

'Ah-v  - =  ° 

i=0  (a  +  c  x) 

sc  ,m  m 


Let  y  =  c  x(m+2) .  Therefore, 
m  m 


z  (j)(-D* 

i=0 


a„  -  y 

l  ,m  m 


a  (N-j+i)  y  nrl-4 

(X  +  — — -  +  -  -SL-) 

U  nH-2  m+2; 


=  0  , 


and  for  large  m. 


J  s 


a  -  y 

v  _ 2 -  =  n 

lo  l  V*/x  Vx 


1=0 
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or 


£  (})(-« 

5=0 


*  * 
xi  ym  xi 1 

1  +  'i(N‘j)  -  -t  +  Hr 


a  C,  /  X 

e  =0 


From  this  equation  we  obtain 
* 


-a  /A  j  x  -a  /A  -a  /A  j-1 


(1  +  -f  (N-j)  -  -r)(l  -  e  “  ) 


•  111  /  HI  \  r\ 

j  e  (1-e  )  =  0  , 


from  which 


or 


y  x.  x_  a  /A 

-f  =  U  +  -y(N-j)  -  -  D1  i 


*  *  a  /A 

y  =  \  +  x  (N-j)  -  x  j/(e  m  -  1)  , 


Vj  =  c  (1+2)  IX  +  Xl(N-j}  "  x’l  j/(e  m  "  1)1  • 


Now  we  note  that  if  j  =0,  then  c  =  —  .  Therefore, 

J  m  a 

m 

X2*,j^[X+Xl  (N-j)  -  x^  j/(eamA-  D] 


for  large  values  of  m. 


(Q.E.D.) 


The  following  table  shows  that  the  approximation  is  close  even  for  small 
values  of  m. 

Incidently,  it  is  readily  seen  that  the  optimal  Bayes  interinspection  time 

X2  j  =  E^9|j’  xl^  is 

J  A  a  *  _m+l 

I  (j)(-l)  [A  +  (N-j+Dx.] 

,  5,-0 

D  - ; — - - - - - - — - - — 

X  j 

2,j  (m-l)  I  cj)(-l>*  [A  +  (N-j+" )x* ] 

5=0  1 


0.5) 


4.  Two  Stage  Optimal  Designs. 
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By  combining  the  results  of  Sections  2  and  3  we  can  construct  two  stage 
optimal  designs  by  using  a  procedure  discussed  by  Zacks  (1973)  and  (1977). 

Now 

1(6, x)  =  N{^|  F(x,0))2{F(x,e)(l-F(x,e))  _1. 

In  order  to  construct  two  stage  optimal  Bayes  design  we  select  n  components, 

0  <  n  <  N  , 

and  perform  our  experiment  at  the  design  level  x,  using  only  n  components. 

Then  at  the  second  stage  the  experiment  is  performed  at  the  design  level  x2  for 
the  remaining  N-n  components.  The  construction  of  the  two  stage  design  involves 
finding  the  vector  (n^,  x^,  x^)  such  that 

(N-n)E{I(0;x2) |F1)  +  nE{I(6;x1)}  £  (N-nQ)E  I(I(e;*J)|F)  +  n0E{I(6;xJ)}  (4.1) 

for  all  n  e  {0,1  ...,N},  x^.x^  e  X. 

That  is  the  use  of  (n^.x^.x^)  gives  global  maximal  Fisher  information  dur¬ 
ing  two  stages. 

We  obtain  (n^.x^.x^)  by  the  following  steps: 

(i)  Find  the  optimal  x^  given  F^  and  n.  Since  x^  is  independent  on  n, 
x^  and  J(x^),  define 

g(n,x1)  =  E[E{I(6;x2)  |  F^}  (N-n)  +n  IOjx^]  . 

(ii)  Determine  n^,  x^  such  that  g(n,x^)  is  maximized. 

(iii)  Redefine  Xq  by  using  n^  and  x^  and  J(x^). 

By  the  lemma  (3.1) 

0  -1  a  A  . 

x2(n,x1> J(Xj))  =  ao)(m+2)  { X  +  (n-J(Xj))xj  -  (e  m  -1)  XjJ(XjH. 


In  order  to  determine  n^,  x^,  one  can  use  the  computer  for  specified  values  of 
N,  A,  m.  If  n  =  0  or  N,  the  solution  is  simple  to  characterize. 

Relative  efficiency  of  the  design  is  defined  as  follows: 

For  the  first  stage,  the  relative  efficiency  RE(o|x,m)  is 

RE  (6  |  A ,  m)  =  I(e;xh/I(0;x*) 

1  m 

where  denotes  the  Bayes  design  level  for  the  first  stage,  i.e., 

x^  =  1.5936  A/m  . 

For  the  two  stage  design  the  relative  efficiency  function  of  0  is  defined  in  a 
similar  manner,  by  using  x^,  where  x^  is  given  by  (3.5). 

For  numerical  values  of  the  relative  efficiency  function  for  different 
values  of  m,  see  Zacks  (1973). 


5.  Kth  Stage  Optimal  Interinspection  Times. 


In  the  present  section  we  consider  the  exponential  failure  distribution. 

If  j, > •  •  •  »j.  denote  the  number  of  components  failing  during  (0,x  ), 

*  *  *  ^  * 

(x^,x^  +  x2),...,(  I  x^,  I  x^),  respectively  then,  we  would  like  to  find  the 

i=1  *  *  * 
(K+l)st  stage  optimal  interinspection  time  x^+^  given  j^,J2* .  •  • ,  jjj*  x^,...,X£. 

.  *  k  k  k 

Now  Jk+1^x' l°(j1(x1) * • • • *Jk^xk> ,xl* # * • ’XjP  B(N,F(x;  ))  where  ^K+1(x)  * 

K  *  * 

number  of  components  failing  during  (  E  x.,  E  x  +  x)  and  the  conditional 

i=l  1  1*1  1 

Fisher  information  of  0  is 


K0;x) 


N  x2  e~x/9 
0*(1  -  e'x/e 


Now  if  0  has  inverted  gamma  (A,m)  as  the  prior  distribution,  then  the  posterior 

k  k  k 

distribution  of  0  given  (jj, j2, . . . , jR,  x1 , x2 , . . . , x^)  Is 
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**■  £  1  X  /  8 

f(e|j)  -  C  E  b(j,;N,F(x,;6))  — jr  e  ;  0  <  6  <  -H=° 

i=i  1  1  em+1 


where  C  is  a  constant.  Therefore, 

~  K 


*  x/fi  m  2  -x/e 

E{I(e,x)  I  j}  =  C  /  I  b(j  ,N,F(x  ;6) )  e  '  ,* . d6 

oh  1  e 


e’(l-e'x/6) 


»  x2  e(M*>/9 


K  „  -  K  -x  /6  j  *  x- 

=  CN  n  (.  )  f  E  [1-e  1  )  1  exp{-(N-i  )x  /6}]  — jt 

i=l  Ji  0  i=l  e 


(1-e 


x/0) 


k  ,  -  K  -x/e  j  . 

=  CN  1  (N  )  x2  /  n  (1-e  ) 


i=l  Ji  0  i=l 
K 


0^5(i-e-x/e) 


exp{-l/6(  E  (N-j.)x.  +  X  +  x) }  d6 . 


i=l 


-x./e  j.  j.  j.  v  * 

Now  (1-e  )  =  E1  (y  )(-l)  exp{-xi  v^e} 

v.=l  i 

l 


Therefore, 


K  -x/0  j. 
fl  (1-e  )  1 

i=l 


K 

n 


3i  1.1 


1=1  v.=l^vi^ 


-  ,  v.  -w  x./e 

E  PU-l)  1  e  1  M 


Ji  ~*2 


JK  i  j 


vrl  v1 


v1 1  11 


ih\ 


K 

E  v. 


k-i) 


i=l 


exp(-  E  v.  x./e) 
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Therefore , 


j  j  j  £ 

E{I(6  ;x)  |  j  }  =  C'  x2  £  7.  ...  7  (  1 )(  2)  (  K )  (-l)i=1 

V1  V2  VK  1  2  K 


-1/0 [Z*v.  x*  +  £*(N-j.)x*  +  X  +  x] 

111  1  Ji  l 


,nrt-5 


x/6. 


de 


0 


0  (1  _  e  '  ) 


K 


9  K  h  K  ji  .  \  Vi 

=  C’  x  (  7  7  )(  I!  (  ))(-l)W 

i=l  v.  i=l  Vi 

i 

x  -1/0  xiI'>1+N-jil  +  X  +  kx] 

£  |  « - — - de 


k=l  0 


0m+5 

■h  ^2  K  j  .  £*v. 


=  C’  x2  £  £  7  ...  XtnC1)} (-1)  1  1 

k=l  v2  v2  vK  i=l  Vi 


1 


(£^x . (v ,+N- j  )  +  X  +  kx)m+^ 

1  l  l  l 


Let 


a  =  x. (v.+N-j , ) 
vi  i'  i  J1 


2 

C"  x  £  £ 

k=l  v . 


JK  K  j . 


.  £  (  £  (,,  )(-!)  X) 


.  ,  v . 
vK  1=1  1 


(£ .  .  a  +  X  +  kx) 
i=l  v . 


Differentiating  with  respect  to  x,  we  get  after  equating  to  zero; 


3x 


E{ I (0 ;x) | j }  =  £  £ 


K  K  j.  v. 

£  (  n  (  X)(-l)  x) 


k=l  u. 


vR  1  Vi 


2(£i_ia  +X)-x(nri-2)k 

i 


(X  +  £^  a  +kx)m<"'* 
1=1  Vi 


=  0 


Let  x  .  be  the  point  for  which  the  above  equation  is  satisfied,  define  a  . 
m,j  m, j 

(x  , ) (m+2) . 
m,  J 


m+4 


Therefore 
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3x 


E{ I (0 ;x) | j }  =  0  impl  ies 


k=l 


jK  K  j  v 

l  [  L  (  )(-l) 

v^=0  i=l  Vi 


2(  A  +  IK  ,  a  )  -  a  .  •  k 

1=1  vt  ro.j 

[(nrt-2)0'.^  a  +.X)  =  ku  .  ]m+"* 

1  =  1  V  IT),  J 


=  0 


a  =x,(v.  +  N  -  )  .  )  . 
v .  l  i  i 

l 


If  all  i  ^  =  0,  1  =  1,2, then  trivially. 


E( I (0 ;x) |  1 }  is  maximized  at  x 


m,  i=0 


where 


m,  i 


i=0 


K  * 
«  (X  +  N  !  .  x.) 

m  _ _  i  =  l 

m  +  2 


We  thus  conjecture  that  if  j  #  0,  then 

jl.  a  (X  +  (N-j .  )x.  - 

*  .01  1=1  'i  i 


,.K 


/ 1 


i-1  Vi/(o 


D) 


m.J 


m  + 


But  this  result  needs  to  be  varified. 


Also,  we  note  that  if  j  =  (0,0,0, ... ,0,  j  ,0,0, ... ,0)  then 

*  ~  am  *  *  an/' 

x  .  ~  — rr  { X  +  (N- j  )x.  -  j .  x.  (e  -1)1 
m,j  nrt-2  i  i  Ji  i 


as  varified  in  Lemma  3.1. 


35 


REFERENCES 

Barlow,  R.E.  and  F.  Proschan  (1967),  Mathematical  Theory  of  Reliability.  John 
Wiley  and  Sons,  Inc.  New  York. 

D'Ortenzio,  R.J.  (1965).  "Approximating  the  Normal  Distribution  Function." 

Systems  Design.  9,  4-7. 

Finney,  D.J.  (1978).  Statistical  Methods  in  Biological  Assays.  3rd  ed. 

Charles  Griffin  and  Corp.  Ltd.  London. 

Freeman,  P.R.  (1970).  "Optimal  Bayesian  Sequential  Estimation  of  the  Median 
Effective  Dose".  Biometrika,  57,  1,  79-89. 

Johnson,  N.L.  and  S.  Kotz  (1970).  Continuous  Univariate  Distributions  I. 

Houghton  Mifflin  Company,  Boston. 

Khan,  M.K.  (1980).  Statistical  Survillance.  Ph.D.  thesis.  Dept,  of  Mathematics 
and  Statistics,  Case  Western  Reserve  llniv.,  Cleveland,  Ohio. 

Widder,  D.  (1961).  Advanced  Calculus.  2nd  ed.  Prentice-Hall,  Inc.  Englewood 
Cliffs,  Nj . 

Zacks,  S.  (1973).  "The  Fisher  Information  Function  and  Design  of  Experiments 

for  Estimation  in  Non-Linear  Statistical  Models".  Technical  Report  No.  8. 
Department  of  Mathemat.es  and  Statistics,  Case  Western  Reserve  L’niv.,  Cleveland. 
Ohio . 

Zacks,  S.  (1977).  "Problems  and  Approaches  in  Design  of  Experiments  for 

Estimation  and  Testing  in  Non-Linear  Models".  Multivariate  Analysis,  P.R. 
Krishnuiah,  ed.  209-223. 

Zacks,  S.  and  W.J.  Fenskc  (1973).  "Sequential  Determination  of  Inspection 

Epochs  for  Reliability  Systems  with  General  Lifetime  Distributions."  Naval 
Res.  Legist.  Quart.  20,  No.  3,  377-385. 


HCuAITy  ClAMIflCATlON  Of  Tull  AAGC  f"»nn  D««  M«iJ) 


REPORT  DOCUMENTATION  PAGE 

RCA V  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  AeClAlSMT'l  CATALOG  NUMICA 

£ _ 

4.  titlc  (m* 

Bayesian  Optimal  Designs  for  Surveillance 
Procedures  in  Attribute  Life  Testing 

1.  rvf  e  Of  ACAOAT  A  ACAIOO  COVCAKO 

Technical  Report 

1.  AUThOAKI 

Mohammad  Kazim  Khan 

1.  CONTAACT  OA  OAAMT  HUMACATfJ 

N000L4-80-C-0325 

10.  AAOGAAA  etCMCNt.  f  AOT«CT.  TA« 

AACA  A  AOAK  UNIT  NUMACAS 

NR  042-276 

M.  CONTAOLLINC  OffICI  NAMC  ano  AOOACfl 

Office  of  Naval  Research 

Statistics  and  Probability  Program  Code  436 
Arlington.  VA  22217  -  - 

*C*OAT  OATt 

December  J5,  1980 

1j.  NUMICR  Of  PAGES 

35 

U.  MONlf OAInG  AGENCY  NAMC  A  AGOAllllll  *//•  Mm  tnm  Conlnllltit  Olllct) 

is.  SECU*ITV  CLASS,  for  t hit  fport) 

mirmrrivMUJi-nu-ru  i  i  1 1 1 1 1  in  mm 

■■■HI 

i*.  oiitaiaution  statcmInt  tot  mi.  n*vtj 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED. 


IT  OIITAIAUTION  1T«TIMII|T  f a I  III*  «A«lracl  MIMlIa  »!•«»  TO.  II  0K»tmt  Ir Mi  "•»•*»> 


1|  iuaauCmCntaay  NOT*! 


l|.  KCY  WOAOI  (Cutlmi*  •"  *14*  II  n»«»«Mrr  Mmrtlr  Or  0l»c*  m»»n) 

Surveillance  Designs,  Fisher  Information,  Inverted  Gamma  Priors, 

Exponential,  logistic,  extreme-value  distributions.  Attribute  Life  Testing 

S( 

'an  |C*MWw*  i  irwiw  <1#*  II  wciiwr  al  lOwllly  Or  nuiA*r? 

Design  of  experiments  for  estimation  of  parameters  in  non-linear  models 
is  studied  in  a  Bayesian  framework,  with  the  objective  of  maximization  of 
the  anticipated  Fisher  information.  Two  stage  optimal  designs  are  proposed 
in  attribute  life  situations. 


00  ,:STn  1473  r.oiTio*  of  '  *ov  »»  ii  oatoiCTt 

:  N  1 1  JJ-  ,c-  Ji  4- 1oG  I  IteuAlTY  C -  »llif  iCATiOm  Of  TH|»  AAGt  Off.  fcwifrfO 


